A Planet in a 0.6-AU Orbit Around the KO Giant HD 102272 



A. Niedzielski 1,2 , K. Gozdziewski 1 , A. Wolszczan 2 , M. Konacki 3,4 , G. Nowak 1 , P. Zielihski 1 



ABSTRACT 



We report the discovery of one or more planet-mass companions to the K0- 
giant HD 102272 with the Hobby-Eberly Telescope. In the absence of any corre- 
lation of the observed periodicities with the standard indicators of stellar activity, 
the observed radial velocity variations are most plausibly explained in terms of 
a Keplerian motion of at least one planet-mass body around the star. With the 
estimated stellar mass of 1.9M , the minimum mass of the confirmed planet is 
5.9Mj. The planet's orbit is characterized by a small but nonzero eccentricity of 
e=0.05 and the semi-major axis of 0.61 AU, which makes it the most compact 
one discovered so far around GK-giants. This detection adds to the existing evi- 
dence that, as predicted by theory, the minimum size of planetary orbits around 
intermediate-mass giants is affected by both planet formation processes and stel- 
lar evolution. The currently available evidence for another planet around HD 
102272 is insufficient to obtain an unambiguous two-orbit solution. 

Subject headings: planetary systems-stars: individual (HD 102272) 



1. Introduction 



Searches for planets around giant stars offer a unique way to extend studies o f planetary 



syste m formation and evolution to stellar masses substantially larger than 1 M Q (ISato et al. 



20031 ). These evolved stars have cool atmospheres and many narrow spectral lines, which can 



be utilized in precision radial velocity (RV) measurements (<10 m s 1 ). Searches for planets 
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aroun d giant stars (jSato et al.ll2008l : lHatzes et al.ll2006l ; iHekker et al.ll2008l ; iNiedzielski et al. 
20071 . and references therein) are beginning to provide the statistics, which are needed to 
constrain the efficiency of planet formatio n as a function o f stel lar m ass and chemical com- 
position. In fact, the initial analyses by lLovis & Mayor (b()07h and iJohnson et al.l (120071 ) 
suggest that the frequency of occurence of massive planets is correlated with stellar mass. Be- 
cause more massive stars probably have more ma ssive disks, these results ap pear to support 
the core accret i on sc enario of planet formation ( Kennedy fc Kenyonl 120081 ). Furthermore, 
Pasquini et al.l (120081 ) have used the apparent lack of correlation between the frequency of 
planets around giants and stellar metallicity to argue that this effect may imply a pollu- 
tion origin of the observ ed planet frequency - metallicity correlation for main sequence stars 
dFischer fc Valentill2005l ). 



In time, the ongoing surveys will also create an experimental basis with wh i ch to 



study the dynamics of p lanetary systems orbiting evolving stars (e.g. iRasio fc Fordl Il996 



Duncan fc Lissauerll 19981 ). Sufficiently large surveys of post-MS giants should furnish enough 
planet detections to meaningfully address the problem of a long-term survival of planetary 
systems around stars that are off the main-sequence (MS) and on their way to the white dwarf 
stage. In fact, the existing data suggest a deficiency or ab sence of orbits with radii below 



0.6-0 .7 AU, which still awaits a fully satisfying explanation (IJohnson et al.ll2007l ; ISato et al. 



20081 ). In addition, the recent detection of a planet around a post-red-giant phase subdwarf 
V391 Peg (jSilvotti et al.ll2007l ) has already demonstrated that planets can survive both the 
RGB and AGB phases of giant evolution. 

In this paper, we describe the detection of a planet around the KO-giant HD 102272 
and discuss a possibility that at least one more planet-mass body exists in this system. An 
outline of the observing procedure and description of the basic properties of HD 102272 are 
given in Section 2, followed by discussion of rotation and stellar activity indicators of the 
star in Section 3. The analysis of radial velocity measurements, is given in Section 4. Finally, 
our results are summarized and further discussed in Section 5. 



2. Observations and properties of the star 



Our surve y, the observing procedure, and data analysis have been described in de- 
tail elsewhere ( INiedzielski et al.l 120071 : INiedzielski fc Wolszczanll2008l ). Briefly, observations 
were made between 2 004, July and 2008, January, with the Hobby-Eberly Telesc ope (HET) 
(IRamsey et al.l Il998l ) equipped with the High Resolu tion Spectrograph (HRS) (lTulllll998l ) 
in the queue scheduled mode (jShetrone et al.l 120071 ) . The spectrograph was used in the 
R=60,000 resolution mode with a gas cell (I2) inserted into the optical path, and it was fed 
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with a 2 arcsec fiber. 



Radial velocities of HD 102272 were measured at 35 epochs spanning the period of 
about 1500 days. Typically, the signal-to-noise ratio per resolution element in the spectra 
was 200-250 at 594 nm in 9-16 minutes of integration, depending on the atmospheric con- 
ditions. Radial ve locities were measured using the standard I2 cell calibration technique 
Butler et al.l (119961 ). A template spectrum was constructed from a high-resolution Fourier 
Transform Spectrometer (FTS) I2 spectrum and a high signal-to-noise stellar spectrum mea- 
sured without the I 2 cell. Doppler shifts were derived from least-squares fits of template 
spectra to stellar spectra with the imprinted I2 absorption lines. The radial velocity for each 
epoch was derived as a mean value of the independent measurements from the 17 usable 
echelle orders with a typic al uncertainty of 6-8 m s _1 at lcr-level. This RV precision made it 
quite sufficient to use the IStumpffi (119801 ) algorithm to refer the measured RVs to the Solar 
System barycenter. 

HD 102272 (BD+14 2434) is a KO-giant with V=8™ 69, B -V=l m 02 and U-B=0 m .69 
( lHaggkvist fc Oja 1 119731 ) and n= 2.76±1.1 1 mas (lESAl Il997l). Its atmospheric param- 
eters were estimated using the method of iTakeda et al.l (J2005aJJb() as T e jj=4908±35 K, 
log(g)=3.07±0.12, and [Fe/Hl=-0 .26±0.08. Compari ng the star's position in the HR di- 
agram with evolutionary tracks of iGirardi et al.l (120001 ). we estimate the mass and radius of 
HD 102272 to be M /M =1.9±Q.3 and R/R =1O.1±4.6, respectively, from calibrations of 
Alonso et all fboooh . 



3. Analysis of the radial velocity data 

Radial velocity measurements of HD 102272 derived from the HRS spectra are shown 
in Figure 1. They reveal a correlated behavior of radial velocity variations of the star on a 
timescale of about 130 days. Moreover, observations made around MJD 54140 indicate that 
the amplitude of the RV curve may vary in time. 



We have used both the nonlinear least-squares and the genetic (GA) algorithm (jCharbonneau 



19951 ) to model the observed RV variations with the standard, six-parameter Keplerian or- 
bits. A fixed 15 m s _1 error was quadratically added to the formal RV uncertainities to 
account for any random RV variations intrinsic to the star, such as fluctuations of the stellar 
surface and possible solar-type oscillations. The anticipated value of this additional jitter 
was conservatively adopted to lie between 20 m s -1 , which is ty pical of stable K-giants 



faekker et al1l2006h . and 5 m s 1 as measured for stable dwarfs (jWrightl 120051 ) . For ex- 



ample the solar-type oscillations alone, when extrapolated for HD 102272 as described by 
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Kjeldsen fc Beddingl (119951 ). would account for a 6 m s" 1 RV variation. The 15 m s" 1 jitter 
is equivalent to f= 0.005 spot on the surface of HD 102272 that would produce 2 m s -1 BVS 
at the same time (IHatzes II2002I ). 



The best-fit of a single orbit and its residuals are shown in Figure 1. The fit, if interpreted 
in terms of an orbiting body, calls for a ~6 Mj up companion in a 127.6-day, slightly eccentric 
orbit around the star (Table 1). The 0.61 A.U. radius of th e orbit is even sm aller than the 
tightest planetary orbit around a red giant identified so far (ISato et al.ll2008l ). However, not 
unexpectedly, the single planet model leaves highly correlated post-fit residuals, indicating 
a possible presence of additional periodicities. 

We have performed a search for the best-fitting two-companion Keplerian model of 
the observed RV variations and identified several solutions of comparable quality with very 
similar x 2 values of the fits. For example, the orbits with the respective approximate periods 
of 179, 350, and 520 days and eccentricities of 0.3, 0.5 and 0.7 are all almost indistinguishable 
in terms of their goodness of fit. At the same time, the parameters of the 127.6-day day 
orbit remain practically unchanged for all of the best two-orbit solutions. An example of 
such a formal solution with the second orbit of 520 days is shown in Figure 1. 

Because of a sparse sampling of the observed RV curve and the resulting ambiguities in 
its modeling, it is useful to perform an additional test of non-randomness of the residuals 
obtained from the best fit of a single planet model to da ta. To accomplish this, we have 
chosen the test of scrambled velocities (IButler et al.ll2004j ) using the GA algorithm. In this 
test, the residuals were randomly scrambled at the epochs of observations and then searched 
for the best-fit parameters of the second orbit. A likelihood, pn, that the residuals represent 
a white noise can be quantified as a ratio of the number of best fits, for which x 2 is l ess 
than or equal to the one derived from the fit to unscrambled data, to the total number of 
trials. The result of the test with 100,000 trials is shown in Figure 2. Clearly, the fit to the 
real signal stands apart from those for the randomly scrambled data with p^ ~ 10~ 5 , which 
strongly suggests the presence of a second, non-random signal in the RV data. 

In principle, the existing ambiguities in the modeling of the RV variations in HD 102272 
can be constrained by imposing the obvious requirement that the true two-planet solution 
is dynamic ally stable. A convenien t way to apply this constraint is to use the GAMP 
algorithm (IGozdziewski et al. 1120061 ). in which an N-body model of the observed RV curve 
has the y 2 criterion mo dified by a penalty term for unstable orbits, and the MEGNO test 
( jCincotta fc Sim6ll2000l ) is applied to check the stability of a given two-planet solution. 



We have performed a GAMP search for the dynamically acceptable orbits of the second 
companion, with the MEGNO integration time set to ~600 orbital periods of the outer 
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companion. The results, in which the x 2 values relative to the best solution found by the 
search are mapped onto the semi-major axis - eccentricity plane, are shown in Figure 3. 

The plot reveals three regions of stability, one of which is located between the mean 
motion resonances (MMR) of 5:2 and 3:1, another one is extended along the 4:1 MMR, and 
yet another one occupies a high-eccentricity region between the 4:1 and 5:1 MMRs. Further 
discussion of these results will be given in the Chapter 5. 



4. Stellar photometry, rotation, and line bisector analysis 



We have followed the standard procedure to verify that the observed RV variations are 
due to orbital motion. It included an examination of the existing photometry of the star for 
possible variations and an analysis of the bisectors and curvatures of the selected spectral 
line profiles. We have also searched for stellar acivity signatures in the Ha line. 

The most precise existing photo metry of HP 102272 consists of 70 photometric obser- 
vations of the star in the Hipparcos JeSA||i9971 ) H p filter between JD 2447877.95490 and 
2448961.44322. The observed scatter in Hp amounted to only 0.015 mag, and no variabil- 
ity above this threshold was detectable. The sta r was also observed by the Northern Sky 
Variability Survey (NSVS) (IWozniak et al.l 120041 ) . We have analyzed the 31 photometric 
observations made between JD=2451318. 188426 and JD=2451630. 300686, which were of a 
sufficient quality, and found no photometric variations of the star. Finally, our Lomb-Scargle 
periodogram analysis of these data did not reveal any spectral peaks above a threshold level 
set at 6 times the mean power computed for periods shorter than one month. Of course, it 
should be kept in mind that the Hipparcos and NSVS data were not simultaneous with our 
RV measurements. 

The projected rotational veloci ty of HP 102272, vs ini =3 ± 1 km s _1 , was estimated 
using the cross-correlation method (IBenz fc Mayor! Il984l ). From this value and the adopted 
stellar radius we have obtained an estimate of the rotation period of P rot ~170 days, which 
is longer than the obs erved 127.6-day period o f radial velocity variation and appears to be 
typical for K gian ts ( de Medeiros et al"1ll996l ). In principle, given the error estimates of 
Alonso et al.l (120001 ). the rotation period of HD 102272 may range from 90 to 250 days. 



Most analyses of the variations of spectral lines using line bisectors are based on a cross- 
correlation function (CCF) representation of an "average" spectral line of the observed star. 
As the I2 lines affect stellar spectra in the iodine cell method of RV determination, they have 
to be properly removed be fore the CCF computation. An efficient method to accomplish 
this has been proposed by (IMartmez Fiorenzano II2005I ). 
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In our implementation of the method, we first cross-correlated the stellar spectrum with 
that of the flat-field I2, in order to measure the wavelength offset between the iodine lines 
present in both spectra. The flat field flux was adjusted to the new wavele ngth scale, adding 
the previously determined offset by using a Hermite spline interpolation (jHill Ill982l ). 



The CCFs were computed by correlating the stellar spectrum with a numerical "zero- 
one", in which the non-zero points were aligned with the positions of stellar absorption 
lines at zero velocity The numerical mask was constructed out of 980 lin es of a very high 



SNR spectrum of HD 17092, the first star with a planet in our survey (INiedzielski et al. 



2007 1), cleaned from spectral features lying within ± 30 km s _1 of the known telluric lines 
Hanuschnik I (120031 ). Out of the 35 resulting bisector measurements, three were affected by 
a low SNR in the original spectra and hence removed from further analysis. The typical 



uncert ainties of the BVS, calculated with the expression given by iMartmez Fiorenzano 
(120051 ) were around ~ 20 m s" 1 . The mean value of BVS was calculated to be 25 ± 21 m 



As both the Ca II K emission line (393.4 nm) and the infrared Ca II triplet lines at 
849.8-854.2 nm are outside the range of our spectra, we used the Ha line (656.28 nm) as a 
chromospheric activity indicator. To minimize a contamination from the telluric lines, we 
measured the EW of the central part of the line profile defined by I/I c <0.7. Our analysis has 
shown that the EW exhibited a 0.69% rms scatter around its mean value of 944.4 ± 6.5mA, 
which was not correlated with the observed RV variability (r = 0.06). 

The analysis of the existing photometry and bisector variability span for HD 102272 
is summarized in Figure 4 which shows the Lomb-Scargle periodograms computed for all 
the time series of interest. We conclude that the only measured parameter, which exhibits 
a significant periodic variability is radial velocity. Similarly, as evident in the plot of BVS 
against the RVs in Figure 5, no correlation between the bisector velocity span and the radial 
velocities exists in HD 102272. 

Using the scatter seen in the Hipparcos photometry of the star, and its rotational velocity 
determined above, we can estimate the amplitude of radial velocity variations an d of the 



bisec tor velocity span due to a possible presence of a spot on the stellar surface (IHatzes 



20021 ). The observed radial velocity amplitude o f HD 102272 is almost an order of magnitude 
larger than 35 m s _1 RV amplitude predicted by lHatzes I (120021 ) for a spot with a filling factor 
of f=0.015, on a star rotating at vsini =3 km s _1 . The expected bisector variations of 5 m 
s _1 are comparable to the precision of our radial velocities measurements and cannot have a 
detectable effect on our results. Similar results have been deri ved with the NSVS phot ometry 
using the spot-induced RV variation modeling described by (jSaar fc Donahue 1119971 ). 
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In principle, it is possible that HD 102272, having no spots over the periods covered by 
both the Hipparcos data (~1083 days) and those from the NSVS (~312 d ays), had devel oped 
a spot later on, at the tim e of our RV measurements. Using the results of lHatzes I (120021 ) and 



Saar fc Donahue I (119971 ). we estimate that, to generate RV variations as large as K=155 m 



s _1 , at a period similar to that of HD 102272 rotation, would require a large spot with the 
filling factor of f=0.17-0.22. Such a spot would translate into a 50-60 m s -1 bisector velocity 
span that is most certainly excluded by our measurements. We also find it unlikely that 
such a large spot, apparently not present in the NSVS data mere ~1400 days prior to our 
observations, would show no evolution in our measurements collected over 1500 days or 11 
cycles of the observed periodicity. Therefore, we conclude that a rotating spot on the surface 
of HD 102272 is unlikely to be a cause of the observed periodic radial velocity variations. 



Discussion 



Our observations of the KO-giant, HD 102272, reveal that the measured radial velocities 
of this star undergo strictly periodic variations at a period of 127.6±0.3 days. When inter- 
preted in terms of a Keplerian motion, this periodicity indicates the presence of a sub-stellar 
companion with a minimum mass of 5.9Mj, in a low eccentricity, e=0.05±0.04 orbit, 0.61 
A.U. away of the star (Table 1). In addition, our data indicate a possible presence of another, 
more distant planetary-mass companion. Unfortunately, in this case, the sparse sampling of 
a putative second orbit does not allow an unambiguous distinction between orbital motion 
and intrinsic stellar effects. In any case, the inner planet has the shortest o rbital period and 



the most compact orbit among planets around GK-giants discovered so far (ISato et al.ll2008 
and references therein). 

As the long period RV variations in red giants may also be related to a combination 
of eff ects including stellar rotation, activity, and non-radial pulsations (e.g. lHatzes et al. 



20061 ). we have analyzed the photo metric data and the behavior of line bisectors of the star, 
following the estab lished practices (Queloz et al.ll200ll ). The details of our procedure have 
been described by ( Niedzielski et al. 20081 ). Its application to HD 102272 has not shown any 
significant correlation between the RV variations and the photometry or line variability of 
the star. Consequently, the most plausible explanation of our data is the presence of at least 
one sub-stellar mass companions around the star. 

The orbital parameters of the confirmed planet are, within errors, independent of the 
choice of the orbit of a putative outer companion. The small orbital radius of the planet 
reinforces the existing evidence that orbits of GK-giant plan ets appear to be wider than ~0.6 
A.U., as suggested by (I Johnson et al.ll2007l ; ISato et al.ll2008l ). These authors have considered 
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the two obvious scenarios to create such a zone of avoidance around the GK-giants, namely 
the tidal capture of a nearby planet by the expanding star, and a paucity of compact orbits 
around giants caused by peculiarities of the disk evolution aroun d intermediate-m ass stars 
( IBurkert &: Ida 120071 ). In particular, simulations carried out by (jSato et al.ll2008l ) indicate 
that, for giants which are over the peak of the red giant branch, a tidal capture of planets 
at orbital radii <0.5 A.U. is possible. Clearly, much more observational evidence is needed 
to place these tentative conclusions on a firmer statistical footing. 

Some of the past reports on planet discoveries around GK-giant s do indicate a p ossible 
existence of more planet-sized companions to the stars in question (jSato et al.ll2008l ) , but 
none of them include analyses designed to constrain their possible orbital parameters. As 
discussed above, and illustrated in Figure 3, our analysis has revealed three regions of plau- 
sible two-planet solutions. The better two of them, located between the MMRs of 5:2 and 
3:1 and close to the 4:1 MMR (Table 1), offer a significant improvement of the quality of the 
fit, in comparison with the one-planet model. In addition, these solutions remain dynami- 
cally stable over at least 1 Gyr, as shown by our numerical integration of the corresponding 
orbits. On the other hand, a more in-depth analysis of these orbits (Gozdziewski et al., in 
preparation) suggests that, because of the compactness of the stability regions for the two 
solutions, a probability for such systems to become trapped in the respective MMRs in the 
course of their dynamical evolution must be very low. Evidently, a confirmation or dismissal 
of a possible dynamical origin of the observed RV variations will have to wait until more 
data become available in the near future. 
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Table 1: Orbital parameters of the two HD 102272 planets derived from the best fit of a 
Keplerian model to the RV data. M is the mean anomaly computed for the epoch of the 
first observation, T = MJD 53042.976890. 



Parameter 


HD 102272 b 


HD102272 c 


P [days] 


127.58 ± 0.30 


520 ± 26 


T [MJD] 


52146 ± 64 


54135 ± 260 


K [ms- 1 ] 


155.5 ± 5.6 


59 ± 11 


e 


0.05 ± 0.04 


0.68 ± 0.06 


u [deg] 


118 ± 58 


320 ± 10 


m,2 sin i [Mj] 


5.9 ± 0.2 


2.6 ± 0.4 


a [AU] 


0.614 ± 0.001 


1.57 ± 0.05 


M(T ) [deg] 


12 ± 60 


324 ± 40 


V [ms- 1 ] 


-94.2 ± 3.6 


(x*) 1/2 


0.87 


a RV [ms 4 ] 


15.4 
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Fig. 1. — (a) Radial velocity measurements of HD 102272 (filled circles) and the best-fit 
models of a single orbit (gray line) and a two-planet system close to the 4:1 MMR (solid 
line), (b-c) The respective best-fit residuals for the above models. 
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Fig. 2. — A histogram of (xl) 1 ^ 2 f° r fits of a single-orbit, Keplerian model to 100,000 sets 
of RV residuals scrambled with a randomly generated signal of the outer planet around HD 
102272. Position of the best-fit solution for the real data is marked with an arrow. 
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Fig. 3. — An ensemble of Newtonian, best-fit models of a two-planet system generated 
with the GAMP code and projected onto the (a c , e c )-plane. A location of each model in 
the diagram is marked with a grey-shaded circle. The consecutive shades of grey (from 
the darkest to the lightest one) correspond to models with the values of (xl) 1 ^ 2 located 
within la, 2a, and 3a of the (xt) 1 ^ 2 °f the best-fit solution identified in the GAMP search, 
respectively. Vertical lines denote positions of low-order mean motion resonances between 
planets b and c. The double-lined curve marks the so called collision line of orbits defined by 
a condition: 05(1 + e&) = a c (l — e c ). As a reference, the best-fit Newtonian model obtained 
without stability constraints is marked with a crossed circle. See text for details. 
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Fig. 4. — Lomb-Scargle periodograms of (a) radial velocities, (b) bisector velocity span, 
(c) Hipparcos H p photometry, and (d) NSVS photometry. The position of the periodicity 
identified in the RV data is indicated with a vertical arrow. Two false alarm probability 
levels, 10% and 1% are marked with dotted lines. 
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Fig. 5. — Radial velocity measurements of HD 102272 plotted against the bisector velocity 
span. Within errors, the BVS remains constant during the period covered by observations 
at the mean value of 25 m s 1 and rms scatter of la = 21 m s _1 . There is no correlation 
between RV and BVS (r=0.03). 



